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ABSTRACT 

We have used the Wisconsin H- Alpha Mapper (WHAM) to observe the spa- 
tially extended distribution of ionized hydrogen in M31 beyond the stellar disk. 
We obtained five sets of observations, centered near the photometric major axis 
of M31, that extend from the center of the galaxy to just off the edge of the 
southwestern HI disk. Beyond the bright stellar disk, but within the HI disk, 
weak Ha is detected with an intensity Iu a = O.OS^q^-R 1 . Since M31 is inclined 
77° with respect to the line of sight, this implies that the ambient intergalactic 
ionizing flux onto each side of M31 is $o — 1-6 x 10 4 photons cm -2 s _1 . Just 
beyond the outer boundary of the HI disk we find no significant detection of 
Ha and place an upper limit lua < 0.01972. 

Subject headings: ISM:general — intergalactic medium — galaxies: ISM — galaxies: 
individual (M31) 



1 1R = Jp- photons cm 2 s 1 ster 1 
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1. Introduction 

The distribution of neutral hydrogen (HI) in spiral galaxies often extends well beyond 
their optical disks (see review by Sancisi & van Albada (1987)). Additionally, some spiral 
systems have neutral hydrogen column densities that sharply truncate at large galactocentric 
distances (van Gorkom 1991; Corbelli & Salpeter 1993). This may be interpreted as the edge 
of the gaseous disk of the galaxy; however, recent theoretical arguments and observational 
evidence suggest that an additional phase of material, ionized hydrogen, may exist beyond 
the regions of observed HI (Maloney 1993; Bland-Hawthorn et al. 1997). Direct observations 
of this ionized gas would yield important physical information about the true size of spiral 
galaxies, and with sufficient sensitivity and spectral resolution, could extend the rotation 
curve of spiral galaxies, yielding more insight into the dynamics of spiral galaxies and their 
distribution of dark matter. Observations of ionized hydrogen beyond the stellar disk may 
also lead to information about the sources of ionization, which may be within the galactic 
system, as in the case of a hot galactic corona, or outside of it. An external source of 
radiation, the intergalactic Lyman continuum flux, plays an important role in cosmological 
models. By observing the brightness of the Ha in the outer HI disk of a local galaxy, one 
may constrain the magnitude of the z = ionizing background flux. 

We have used the Wisconsin H-Alpha Mapper (WHAM) to observe the distribution of 
diffuse ionized hydrogen in the local spiral galaxy M31. The high sensitivity and large field 
of view of WHAM make M31 the best galaxy in which to search for an extended disk of 
ionized hydrogen. 

2. Observations 

The WHAM instrument is a fully remotely operated facility with a 15 cm, dual-etalon 
Fabry-Perot spectrometer at the focal plane of a 0.6 m telescope atop Kitt Peak in Arizona 
(Reynolds et al. 1998; Tufte 1997). The WHAM spectrometer has a 1° diameter circular field 
of view on the sky, and a velocity resolution of 12 kms -1 within a 200 kms -1 wide spectral 
window that can be centered on any wavelength between 4800 A and 7300 A. WHAM was 
designed to detect very weak emission lines from ionized gas. 

For the observations of M31, the spectral window was centered at approximately -500 
kms -1 with respect to rest Ha, which is the approximate radial velocity of the southwestern 
half of the galaxy. Spectra were obtained toward a total of nine lines of sight. Five lines of 
sight ("ONs") were centered near the photometric major axis of M31 and extended from the 
center of M31 to just off the edge of the observed southwestern HI disk (see Figure 1). Four 



- 3- 



lines of sight ("OFFs"), located a couple of degrees away from M31 (see Table 1) and free 
of M31 Ha emission, were used as comparison spectra. Each direction was observed for 300 
s at a time, alternating between OFFs and ONs. The data were obtained over a total of five 
nights during three different periods separated by almost one year. This was done to take 
advantage of the motion of the expected M31 Ha emission line with respect to fixed weak 
atmospheric lines. Table 1 summarizes the details of the observations. 

The WHAM interferometer produced 128x128 pixel ring images from which spectra 
were extracted via an annular ring summing technique; each position within the 1° beam 
is equally sampled by each spectral element (Tufte 1997). The resultant spectra were flat- 
fielded, continuum subtracted, and intensity calibrated. Signatures of cosmic rays and phys- 
ical defects in the interference filter were removed with pixel masking. 

Figure 2 shows an average of 44 spectra toward all OFF lines of sight taken on the same 
night, free of M31 emission. The bright feature in the red side of the spectrum is a well- 
known OH atmospheric line at 6553.6 A, which we used for precise wavelength calibration. 
The feature near -480 kms -1 is a Fabry-Perot ghost from the bright geocoronal Ha line near 
kms^ 1 ; the other features are extremely weak, unidentified atmospheric lines within the 
spectral window. The strengths of all of these lines vary both temporally and spatially on 
the sky, which make the task of identifying the weak Ha emission line from M31 formidable. 

To isolate the emission solely due to M31, the spectra obtained in the OFF directions 
were subtracted from the ON spectra. Due to the complex nature of the intensity variations 
of the weak atmospheric features, only a subset of the OFF spectra were subtracted from 
the ONs. This subset satisfied the criteria that the resultant spectra formed a flat baseline 
around the region of expected M31 emission. This process produces accurate spectra if 
the OFFs possess no significant extra-terrestrial emission features within the spectrum. The 
potentially incomplete subtraction of the atmospheric lines is the largest source of uncertainty 
in our results. 



3. Results 

Beam 1, positioned at the center of M31, reveals significant Ha emission. However, 
due to the dynamical location of the line of sight, the width of the line extended beyond the 
±100 kms -1 range of the spectrum. Hence, an accurate spectral baseline and estimate of the 
emission line strength was not possible. For the remaining ON beam positions, the expected 
emission line width should be relatively narrow, comparable to the thermal and turbulent 
motion of the emitting gas. This is confirmed by observations of M31's flat rotation curve 
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toward these directions, as seen in the 21 cm spectrum along the major axis (Cram et al. 
1980). 

The spectra taken toward Beam 2 also reveals significant Ha emission, as seen in Figure 
3, with an intensity of l Ha = 0.66^/?. The error in this measurement is due to the 
uncertainty in the spectral baseline generated by the incomplete subtraction of the weak 
atmospheric features. Toward the line of sight off of the stellar disk, but within the HI disk 
(Beam 3), Ha emission is detected near the velocity of the HI, with l^a — 0.051qo2-^ ( see 
Figure 4). Figure 1 shows that this emission arises far from the optical stellar disk. 

Toward the two lines of sight just off the edge of the HI disk (Beams 4 & 5), no significant 
Ha emission was detected. Upper limits of In a < 0.019i? and In a < 0.031i? were determined 
for Beams 4 and 5, respectively (see Figure 4). These are the strengths of the strongest 
emission that can be placed at the HI velocity, with the same width as the Ha emission 
detected in Beam 3, and not be discernable in the spectra. Note that the scatter in the 
spectra toward Beam 5 is considerably larger than for Beam 4. This increase in the scatter 
arises from intrinsically larger variations in the atmospheric lines in this direction. A more 
accurate removal of these lines toward Beam 4 has provided a stronger upper limit, which 
we adopt as the largest Ha emission strength allowable off the edge of the HI disk in M31. 
Although the data indicate a somewhat lower limit, we have concluded that the variations 
in the atmospheric lines, which produce incomplete subtraction of these features in the ON- 
OFF spectra, could hide a signal up as large as 0.019.R. We note that the 2a uncertainty 
due to Poisson noise in the data is much lower, at 0.004/?. 

Other deep CCD imaging searches for faint Ha emission in external galaxies (Zurita 
et al. 2000; Ferguson et al. 1998a; Martin & Kennicutt 2001; Walterbos & Braun 1992) 
are sensitive to signals ranging from 0.5 to 100i?, with typical sensitivities of tens of 
Rayleighs. Our ability to detect emission down to ~ O.Oli? is a result of the unique dual- 
etalon, high resolution Fabry- Perot design of WHAM, combined with the power of the ON- 
OFF technique. 

4. Discussion and Conclusions 

One of the original motivations for these deep Ha observations of M31 was to search for 
HII extending beyond the radius of the HI cutoff. No such extended emission was found down 
to our detection limit of 0.019/2. This corresponds to a face-on emission measure EM± = 
0.016 cm -6 pc for a temperature of 10 4 K, and a column recombination rate of 1.2 x 10 4 
hydrogen recombinations cm -2 s -1 . To calculate the emission measure and recombination 



- 5 - 



rate we estimated the extinction of the Ha photons as they pass through the Milky Way 
along each of the lines of sight. We used the Galactic neutral hydrogen column density Nhi 
toward each direction (Hartmann & Burton 1997), and the relationship between Nhi and t\ 
given by Bohlin et al. (1978) and Mathis (1990) to increase ln a by ps 30% over the observed 
ln a - The extinction within the outer disk of M31 is negligible as compared to the foreground 
Galaxy (Cuillandre et al. 2001), and we therefore have not corrected for that effect. 

On the other hand, we did detect ionized hydrogen in the extended HI disk of M31 
(Beam 3). One possible source of this ionization is young, massive O stars within the HI 
layer. Cuillandre et al. (2001) recently completed a deep, V and / band photometric imaging 
study of a region within our Beam 3, and did find evidence of young B stars. However, they 
also mention that while an estimate of the threshold for very massive star formation in 
this region is difficult to ascertain, it is unlikely that the gas disk surface density exceeds 
the critical density (Martin & Kennicutt 2001). Therefore, while deep Ha imaging studies 
of galaxies of later-type than M31 (Sb) have revealed sites of massive star formation at 
comparably large galactocentric distances (Ferguson et al. 1998a,b), it seems unlikely that 
the observed Ha emission toward Beam 3 arises from photoionization from massive O stars. 

Another possible source of this HII is an external flux of Lyman continuum radiation. 
Since the HI fills the WHAM beam for position 3 (see Cram et al. (1980)), the observed 
Ha intensity of 0.05i2 provides a direct estimate of the required ionizing flux (Maloney 1993). 
Assuming that the HI absorbs all of the incident Lyman continuum photons and that the gas 
is optically thin to Ha, we find that an ionizing flux incident $o — 1-6 x 10 4 photons cm -2 s _1 
on each side of M31 would account for the Ha. We derived this value assuming case B and 
an electron temperature T e = 10 4 K, making the small extinction correction described above, 
and taking into account the 77° inclination of M31's disk from the line of sight. However, 
since additional sources of ionization may be present, e.g. B stars, planetary nebulae, or hot 
white dwarf stars in the outer disk and halo of M31 (Cuillandre et al. 2001), this value of $ 
must be considered an upper limit on the external ionizing flux. Other authors have derived 
$o with various techniques combining both observations and theory and probing a variety 
of extragalactic environments (see Table 1 of Shull et al. (1999)). We find that our value of 
$o = 1.6 x 10 4 photons cm -2 s -1 is among the lowest of these reported values. 

Our upper limit on $o near M31 implies that a spherical, isolated HI cloud bathed only 
in an isotropic intergalactic radiation field of the Local Group would have an Ha intensity 
Ihq. < 0.03/2. This is significantly (~ factor of 3) below the observed intensities of high 
velocity HI clouds in the Milky Way's halo (Tufte et al. 1998) and in the Magellanic Stream 
(Weiner & Williams 1996), implying an additional source of ionizing radiation in the halo of 
the Milky Way, such as Lyman continuum photons from O stars escaping the disk (Bland- 



- 6- 



Hawthorn & Maloney 1999; Dove, Shull, & Ferrara 2000). Ha observations of high velocity 
HI clouds may therefore provide important information about the location of clouds within 
the Local Group. Specifically, clouds far from a galaxy (e.g. Blitz et al. (1999)) should have 
Ha intensities < 0.03-R. 

GJM, RJR, LMH, and SLT acknowledge support from the National Science Foundation 
through grant AST 96-19424. PRM is supported by the National Science Foundation under 
grant AST 99-00871. 
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Table I. Summary of Observations 
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a Coordinates given in [B1950] equinox. 

b Galactocentric distance, using a distance of 784 kpc to M31 (Stanek & Garnavich 1998). 
c Radius of the fiB = 25 mag arcsec~ 2 isophote (de Vaucouleurs et al. 1991). 
d Velocity of the Ha emission with respect to the LSR. 



Fig. 1. — Overlay of an optical image of M31, HI contours, and the ON lines of sight 
considered in this study. The optical image is from the DSS-II at STScI. The HI contour 
diagram is adapted from Cram et al. (1980). The outermost contour level is at N#j = 
6.51 x 10 19 cm~ 2 . Detectable HI emission extends beyond this last contour to fill Beam 3. 
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Fig. 2. — Average spectrum of 44 OFF lines of sight, all taken on the same night, for a total 
integration of 3 h, 40 min. The y-axis is given in milli-Rayleighs per kms -1 . The top graph 
displays the total spectrum; the bottom graph is an expansion of the full spectrum at low 
intensity levels. Note the strong OH line near -410 kms -1 , the geocoronal Ha ghost at -480 
kms -1 , and several weaker atmospheric lines. 
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Fig. 3. — Reduced spectrum of the region toward the southwestern stellar disk (Beam 2). 
The M31 Ha emission at -535 kms -1 has a strength of ln a = 0.66+o;ogi?. The data are 
represented as diamonds and the best fit by a solid line through the data. Note the incomplete 
subtraction of the OH line, and the location of the continuum. The straight, solid line is the 
continuum used for the best fit, whereas the dashed and dotted lines are upper and lower 
limits to the continuum, reflected in the upper and lower error limits of ln a - 
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Fig. 4. — Reduced Ha spectra away from M31's stellar disk toward Beams 3, 4 and 5. 
Ha emission from M31 is present in the top plot toward Beam 3 at a velocity of -521 kms -1 , 
with a strength of l Ha = O.OS+q^-R (symbols and lines as in Fig. 3). The middle and bottom 
plots of spectra towards Beams 4 & 5 show no detectable Ha emission. The solid line in 
these spectra are reasonable lower limits to their continua plus an Ha emission line with the 
same velocity and width as the emission toward Beam 3. Lines of strength = 0.019/? 
and Ina — 0.031i? are overlayed in Beams 4 and 5, respectively. Note the scatter in the data 
and that the vertical scales in all three plots are identical. 



